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SECTION  1 


INTRODUCTION 


As  a  result  of  several  recent  programs  over  the  past  few  years,  MITRE  developed 
software  to  predict  the  performance  of  airborne  phased  array  radars  and  algorithms  using 
space-time  adaptive  processing  (STAR)  to  suppress  interference  signals  received  by  these 
radars.  There  are  two  types  of  software  capabilities,  one  is  the  steady-state  performance 
prediction  simulation,  which  can  model  the  environment  of  interference  signals  and  the  other 
is  finite  sample  processing,  which  can  process  measured  or  simulated  in-phase  and 
quadrature  (I  &  Q)  data.  These  software  capabilities  have  not  been  adequately  evaluated 
with  measured  data.  Rome  Laboratory’s  (RL’s)  Multichannel  Airborne  Radar  Measurements 
(MCARM)  program  is  currently  developing  a  testbed  that  will  obtain  this  type  of  data.  The 
objective  of  this  project  is  to  enhance  MITRE’s  simulation  and  processing  capability  to 
develop  benchmarks  for  the  RL  testbed  for  determining  its  potential  performance,  and  to 
process  and  evaluate  the  measured  data  that  will  be  obtained  from  this  testbed  in  FY95. 

The  project  can  be  described  by  four  overall  tasks.  Two  tasks  cover  upgrades  to  each  of 
the  software  capabilities.  One  is  the  application  of  an  improved  graphics  interface,  and  one  is 
a  demonstration  of  the  steady-state  simulation  plus  recommendations  for  their  application. 

•  Several  enhancements  have  been  made  to  the  capability  of  the  simulation  model  that 
predicts  the  steady-state  performance.  Elevation  degrees  of  freedom  (DOF)  and  a 
model  for  multiple  point  near-field  scatterers  has  been  added.  Also,  the  outputs  of  a 
simulation  code  based  upon  geometric  optics  and  diffraction  that  provides  more 
detailed  models  of  the  interaction  with  aircraft  scatterers  can  be  in.serted  directly  into 
the  simulation.  The  simulation  has  been  reorganized  to  run  more  efficiently. 

•  The  finite  sample  processing  software  can  now  apply  spatial  and/or  temporal 
preprocessing  transformations  to  a  cube  of  received  I  &  Q  data.  The  covariance  matrix 
used  to  calculate  the  adapted  weights  is  obtained  from  the  transformed  data. 

•  A  mouse  driven  front-end  with  menus  has  been  applied  to  the  simulation  and 
processing  software.  A  graphical  user  interface  (GUI)  with  several  plotting  package 
options  has  been  provided  for  showing  the  input  files  and  the  output  performance 
measures.  Most  of  these  have  cross-platform  capability. 

•  The  simulation  has  been  demonstrated  for  .scenarios  based  upon  the  planned  MCARM 
experiments.  Future  applications  of  these  simulations  have  been  recommended  based 
upon  these  experiments,  potential  Airborne  Warning  and  Control  System  (AW ACS) 
applications,  the  ARPA  (Advanced  Research  Projects  Agency)  mountain  top  data 
reduction,  and  the  proposed  Bistatic  Adjunct  Surveillance  System  (BASS). 
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SECTION  2 


STEADY-STATE  SIMULATION 


mitre’s  Advanced  Airborne  Radar  Simulation  is  based  upon  the  STAR  architecture  shown 
in  Figure  1  [1-4].  Tapped  delay  lines  are  placed  at  the  output  of  each  array  output,  with  the 
taps  spaced  by  one  pulse  repetition  interval  (PRI).  The  adapted  weights  are  obtained  by 
multiplying  the  inverse  of  the  covariance  matrix  created  from  the  received  signals 
(interference-plus-target-plus-noise  signals)  at  each  of  the  tap  outputs  with  a  steering  vector. 
Each  weight  is  applied  to  the  appropriate  PRI  tap  (i.e.,  pulse),  and  the  weighted  pulse  outputs 
are  coherently  summed  to  provide  the  adapted  output  signal.  The  number  of  DOF  is  the 
number  of  pulses  times  the  number  of  array  outputs  (elements,  beams,  or  subarrays).  When 
the  number  of  pulses  used  in  the  STAP  filter  is  less  than  the  number  of  PRI  outputs  to  be 
coherently  processed,  Doppler  processing  following  the  STAP  filter  summation  provides 
further  coherent  gain  for  the  target  signal.  The  STAP  performance  (e.g.,  signal-to- 
interference  rado  (SIR))  can  be  measured  at  the  output  of  the  summadon  and  at  any  Doppler 
filter  output.  Additional  temporal  taps  spaced  by  a  sampling  interval  that  is  less  than  the 
length  of  a  compressed  pulse  represents  another  dimension  of  DOF  that  can  be  applied  to  the 
output  of  each  PRI  tap  for  suppressing  multipath  and  for  channel  mismatch  [4].  The 
individual  steady-state  covariance  matrices  for  the  received  clutter  signals,  jammer  signals, 
and  thermal  noise  are  calculated  based  upon  their  modeled  spatial  and  temporal  correlation 
properties  as  affected  by  the  platform  velocity  and  crab  angle,  internal  clutter  motion,  signal 
bandwidth,  and  the  match  of  the  output  channels.  Additional  insight  into  the  adapted 
performance  is  provided  by  calculating  the  eigenspectra  of  the  interference  covariance  matrix 
and  the  adapted  antenna  gain  patterns.  Figure  2  is  a  global  flow  chart  showing  the  options  of 
the  steady-state  simulation  with  the  shaded  portions  identifying  the  new  capabilities  that  have 
been  provided. 


2.1  ELEVATION  DOF 

The  Advanced  Airborne  Radar  Simulation  employs  a  planar  array  with  a  specified  aperture 
weighting  and  with  the  antenna  element  gain  pattern.  Previously,  the  simulation 
characterized  each  column  of  the  planar  array  on  receive  as  a  single  non-adapted  elevation 
subarray.  The  effects  of  azimuth  spatial  DOF  on  the  performance  of  STAP  can  be  evaluated 
by  adaptively  proce.ssing  the  outputs  of  the.se  receive  azimuth  subarrays.  However,  for 
medium-  and  high-pulse  repetition  frequency  (PRF)  waveforms,  the  long-range  target  and 
short-range  clutter  may  compete  with  each  other  becau.se  of  range  folding.  For  the.se 
waveforms,  the  space-time  processing  (STP)  architecture  may  require  additional  spatial  DOF 
in  elevation  to  reduce  the  contributions  of  the  .short-range  clutter  as  discu.s.sed  below. 
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Element  1  Element  2  Element  E 


Doppler  Filter  Outputs 


Figure  1.  Jointly- Adaptive  Space-Time  Processing  Architectures 
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The  simulation  has  been  modified  to  allow  the  caleulation  an  equal  number  of  multiple  non- 
adapted  subarrays  in  elevation  for  each  array  column  or  azimuth  subarray.  Earlier  the 
simulation  had  52  subroutines.  Fifteen  subroutines  were  modified  to  incorporate  this 
adaptive  elevation  DOF  capability. 

Two  examples  illustrate  the  elevation  DOF  feature  of  the  simulation.  First,  we  validate  the 
program  using  a  horizontal  linear  an-ay  with  a  .single  sidelobe  Jammer  in  a  clutter-free 
environment.  We  then  rotated  the  array  by  90  degrees  and  placed  the  single  jammer  at  the 
same  angular  location  in  elevation.  The  jammer-plu.s-noise-to-noise  ratio  (JNR)  results  from 
the  two  programs  were  identical  when  using  the  modified  simulation.  Second,  we  illustrate 
how  elevation  DOF  can  be  u.sed  to  mitigate  short-  and  long-range  clutter  competing  with  a 
long-range  target.  Although  the  spatial  elevation  DOF  alone  can  mitigate  short-range  clutter, 
additional  temporal  DOF  are  required  to  cancel  long-range  clutter  while  maintaining-  the 
mainbeam  and  not  canceling  the  target. 

2.L1  Sidelobe  Jammer  Nulling  Example 

In  this  example,  we  chose  a  horizontal  20-element  linear  array  electronically  .scanned  in 
azimuth  to  -10  degrees  (relative  to  broadside)  and  a  sidelobe  jammer  located  in  azimuth  at 
-32  degrees  (relative  to  broadside).  Using  the  previous  steady-state  version  of  the  simulation, 
the  quiescent  JNR  was  44  dB,  and  the  adapted  JNR  using  element  space  processing  (i.e., 

20  azimuth  DOF)  was  8.2  dB.  Note  that  the  jammer  was  not  completely  canceled  due  to  a 
specified  cancellation  ratio  of  55  dB.  Next,  a  vertical  20-element  linear  array  (i.e.,  20 
elevation  DOF)  electronically  scanned  in  elevation  to  10  degrees  was  used  to  cancel  a 
sidelobe  jammer  at  32  degrees  in  elevation.  Figure  3  shows  the  quiescent  and  adapted  JNR 
using  the  modified  simulation.  Figure  4  shows  the  quiescent  and  adapted  elevation  receive 
antenna  gain  patterns,  illustrating  that  there  is  an  adapted  null  at  32  degrees  used  to  cancel  the 
sidelobe  jammer.  The  performance  mea.sures  and  patterns  were  identical  for  the  horizontal 
and  vertical  aixay  configurations  demonstrating  the  successful  implementation  of  two 
dimensions  of  spatial  DOF  in  the  simulation. 

2.1.2  Clutter  Suppression  Example 

This  example  illustrates  how  elevation  DOF  can  adaptively  cancel  short-  and  long-range 
clutter  competing  with  a  long-range  target.  For  this  example  we  use  a  30  x  8  planar  array 
electronically  .scanned  in  elevation  to  2.8  degrees.  The  PRF  of  the  radar  waveform  is  1400 
Hz,  which  cau.ses  an  additional  three  range  ambiguities  at  elevation  angles  of  3.1, 4.9,  and 
75.4  degrees,  respectively.  Figure  5  shows  the  performance  measures  when  eight  elevation 
DOF  are  u.sed  to  adaptively  cancel  the  clutter.  The  dutler-plus-noi.se-to-noi.se  ratio  (CNR)  is 
reduced  from  60  to  0.3  dB.  However,  there  is  also  target  cancellation,  as  indicated  by  the 
significant  loss  in  .signal-to-noi.se  ratio  (SNR)  and  small  improvement  in  SIR.  The  quie.scent 
and  adapted  elevation  receive  antenna  gain  patterns  shows  the  cancellation  of  both  the  target 
and  clutter.  Figure  6  illustrates  that  an  adaptive  null  in  the  elevation  gain  pattern  at  75.4 
degrees  cancels  the  short-range  clutter. 
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Measures  (in  dB)  Adaptive  Quiescent 


J  +  N 
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37.09 
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J  +  N 
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=  jammer-plus- noise- to- noise  ratio 


=  signal-to-interference  ratio 


=  signal-to-noise  ratio 


44.02 


-6.91 


37.11 


Figure  3.  Performance  Measures  for  Canceling  a  Sidelobe  Jammer  With  20 

Elevation  DOF 
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Number  of  Antenna  Columns  =  1  Radar  Bandwidth  =  1 .0E06  Hz 

Number  of  Antenna  Rows  =  20  Jammer  Model  =  Wideband 
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Figure  4.  Elevation  Receive  Antenna  Gain  Patterns  for  Canceling  a  Sidelobe 

Jammer  with  211  Elevation  DOF 
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Performance  Measures  for  Clutter  Cancellation  with  8  Elevation 

DOF 
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Figure  6.  Elevation  Receive  Antenna  Gain  Patterns  for  Clutter  Cancellation 

with  8  FJIevation  DOF 


However,  there  are  adaptive  nulls  in  the  mainlobe  corresponding  to  the  elevation  angle 
locations  of  the  long-range  clutter  (which  includes  the  elevation  location  of  the  target). 

Temporal  DOF  (i.e.,  pulses)  can  adaptively  discriminate  a  moving  target  from  the  long-range 
clutter  in  the  main  beam.  Figure  7  shows  the  performance  measures  using  32  DOF  where 
four  temporal  and  eight  elevation  DOF  are  used  simultaneously  to  cancel  the  clutter. 
Eigenvalue  compensation  was  applied  to  minimize  antenna  pattern  distortion  [3].  The  CNR 
is  reduced  from  60.0  to  6.5  dB,  while  the  SNR  is  increa.sed  from  5.7  to  10.7  dB,  resulting  in  a 
significant  improvement  in  SIR  from  -54.4  to  4.1  dB.  Figure  8  illustrates  the  elevation 
receive  antenna  gain  pattern  evaluated  at  the  target  Doppler  (i.e.,  half  the  blind  speed). 


2.2  NEAR-FIELD  SCATTERING 

Near- field  scattering  by  the  radar’s  platform  can  degrade  the  free-space  pattern  of  the  antenna 
by  redirecting  energy  from  the  sidelobes  into  the  mainbeam,  thereby  degrading  the  radars 
performance.  STP  has  been  shown  to  be  an  effective  technique  for  mitigating  near-field 
scattering  [6].  For  each  clutter  scatterer,  the  software  calculates  an  entry  of  the  steady-state 
covariance  matrix  based  on  the  radar  range  equation,  the  free-space  transmit  gain  pattern  of 
the  array,  and  the  linear  phase  terms  due  to  the  clutter’s  spatial  (i.e.,  element-to-element)  and 
temporal  (i.e.,  tap-to-tap  and  pulse-to-pulse)  correlation  properties.  In  the  absence  of  near- 
field  scattering,  the  clutter  steady:State  covariance  matrix  is  spatially  and  temporally 
stationary,  and  is  constructed  by  calculating  only  a  single  row  of  entries  and  exploiting  its 
Toeplitz  structure.  However,  one  of  the  critical  factors  that  can  limit  the  performance  of 
STAP  is  the  near-field  scattering  effects  due  to  the  antenna-aircraft  interactions.  The 
modifications  to  include  the  scattering  effects  from  near-field  obstacles  in  the  simulation’s 
steady-state  mode  are  described  below. 

2.2.1  Point  Scatterer  Model 

The  effects  of  multipath  caused  by  point  scatterers  in  the  near-field  of  the  array  antenna  can 
be  described  by  the  sum  of  four  .scattering  paths  as  illustrated  in  Figure  9.  The  direct  transmit 
and  receive  path  between  the  airay  and  a  clutter  patch  in  the  far  field  of  the  array  on  the 
ground  is  shown  in  Figure  9a.  Near-field  .scattering  creates  three  additional  paths.  One 
bounce  path  on  receive  as  in  Figure  9b  requires  augmenting  the  steady-state  covariance 
matrix  calculation  with  additional  non-linear  pha.se  ternis  due  the  nonstationary  spatial 
correlation  properties  of  the  received  clutter.  One  bounce  path  on  transmit  as  in  Figure  9c 
requires  the  calculation  of  the  far-field  transmit  gain  pattern  in  the  pre.sence  of  the  near-field 
point  .scatterers.  The  two-bounce  path  in  Figure  9d  is  a  composite  of  one  bounce  path  on 
transmit  and  one  on  receive.  These  three  added  bounce  paths  cause  the  steady-state 
covariance  matrix  for  the  received  clutter  to  be  non-Toeplitz,  where  all  entries  of  the  matrix 
need  to  be  calculated  explicitly.  Therefore,  this  modification  can  significantly  increa.se  the 
computational  complexity  and  run  time  of  the  simulation. 
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Figure  7.  Performance  Measures  for  Clutter  Cancellation  with 
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Figure  9.  Four  Simulated  Scattering  Paths 


2.2.2  Examples  of  Point  Scatterer  Model 

Two  examples  are  presented  to  illustrate  the  near-field  point  scattering  model  of  the 
simulation.  First  is  the  effect  of  near-field  scattering  on  the  cancellation  of  a  single  .sidelobe 
jammer  and  second  is  its  effect  on  clutter  cancellation. 

A  16-element  horizontal  linear  array  was  simulated,  with  the  antenna-platform  interactions 
represented  by  a  line  of  eleven  closely  spaced  point  .scatterers  (each  with  a  bistatic  cross 
section  of  2  m^),  as  shown  in  Figure  10.  The  array  is  collinear  with  the  x-axis  and  centered 
about  the  origin,  and  the  platform  heading  is  in  the  positive  x-direction.  The  single  sidelobe 
jammer  is  located  at  -43  degrees  in  azimuth.  Figure  1 1  shows  the  quiescent  and  adapted 
receive  antenna  gain  patterns  without  near-field  scattering.  Figure  1 1  shows  the  quie.scent 
and  adapted  performance  measures  in  the  absence  of  near-field  scattering.  The  jammer  was 
not  completely  canceled  (i.e.,  adapted  JNR  equal  2.9  dB)  due  to  the  channel  matching, 
cancellation  ratio  of  50  dB  as  shown  in  Figure  12.  In  comparison.  Figure  13  shows  the 
quiescent  and  adapted  patterns  in  the  presence  of  the  near-field  point  .scatterers.  The 
quiescent  gain  in  the  azimuth  direction  of  the  jammer  increased  the  JNR  by  almost  13  dB,  as 
show  in  Figure  14.  The  near-field  .scattering  and  channel  mismatch  effects  increa.sed  the 
adapted  JNR  to  6.3  dB.  However,  if  the  bistatic  cross  .section  of  each  of  the  near-field 
scatterers  is  reduced,  the  adapted  performance  approaches  the  level  achieved  without  near¬ 
field  scattering. 

The  effect  of  the  line  of  point  near-field  scatterers  on  clutter  suppression  was  illustrated 
with  a  horizontal  linear  array  of  columns.  Figure  15  shows  the  adapted  SIR  without  near¬ 
field  scattering  for  a  varying  number  of  spatial  DOF  or  columns  and  temporal  DOF  or  pulses, 
The  adapted  SIR  is  normalized  by  the  adapted  SIR  obtained  in  the  noise-only  case  resulting 
in  a  maximum  achievable  value  of  0  dB.  From  Figure  15,  for  example,  canceling  the  clutter 
with  an  internal  motion  of  0.1  m/s  and  adapted  SIR  of -3  dB  requires  eight  elements  and 
three  pulses.  Figure  16  shows  the  adapted  SIR  for  the  same  array  with  the  line  of  eleven 
near-field  point  scatters.  In  order  to  achieve  an  adapted  SIR  of  less  than  3  dB  with  eight  to 
twelve  columns  now  requires  four  pulses  and  there  is  a  considerable  degradation  with  only 
four  columns. 

2.2.3  Geometric  Theory  of  Diffraction  (GTD)  Model 

The  .second  modification  uses  a  simulation  that  models  the  electromagnetic  interaction 
between  an  antenna  and  its  platform  (e.g.,  the  airplane  on  which  it  is  mounted).  The 
Numerical  Electromagnetics  Code-Basic  Scattering  Code  (NECBSC),  a  simulation  that  is 
ba.sed  upon  geometric  optics  and  the  GTD.  It  deieiTnines  the  complex  amplitude  and  pha.se 
information  at  each  far-field  point  (corresponding  to  a  clutter  .scatterer  location)  due  to  each 
transmit  antenna  element  in  the  pre.sence  of  a  model  of  the  aircraft  structure.  This  .simulation 
was  developed  by  Ohio  State  University  to  model  interactions  between  antennas  and  their 
platforms  [7].  The  signal  at  each  of  the  array  outputs  obtained  from  this  GTD-ba.sed 
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Figure  10.  Location  of  Eleven  Near-Field  Point  Scatterers 
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Figure  11.  Azimuth  Receive  Antenna  Gain  Patterns  for  Cancelling  a  Sidelobe  Jammer,  Without  Near-Field 

Scattering 
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Figure  12.  Performance  Measures  for  Cancelling  a  Sidelobe  Jammer  Without  Near- 

Field  Scattering 
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Number  of  Antenna  Columns  =  1 6  Jammer  Model  =  Wideband 

Number  of  Antenna  Rows  =  1  Jammer  Azimuth  Angle  =  -  43  Degrees 
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Figure  13.  Azimuth  Receive  Antenna  Gain  Patterns  for  Cancelling  a  Sidelobe  Jammer,  With  Near-Field 

Scattering 


182472 


Performance  * 

Measures  (in  dB)  Adaptive 


(LLE) 

.  N  , 


6.35 


16.97 


23.31 


J  +  NA 
N  . 


=  jammer-plus-noise-to-noise  ratio 


=  signal-to-interference  ratio 


=  signal- to- noise  ratio 


Quiescent 

51.61 


-27.09 


24.52 


Figure  14.  Performance  Measures  for  Cancelling  a  Sidelobe  Jammer  With  Near- 

Field  Scattering 
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Figure  15.  Performance  Measures  for  Cancelling  Clutter,  Without  Near-Field  Scattering 
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simulation  is  stored  on  an  output  file  that  can  be  read  into  the  steady  state  STAP  simulation 
to  construct  the  clutter  steady-state  covariance  matrix  [8].  The  STAP  performance  predicted 
using  this  modification  is  based  upon  a  more  accurate  model  of  both  the  near-field  scatterer 
and  the  electromagnetic  scattering  phenomenology  in  comparison  with  the  point  scatterer 
model.  However,  this  is  provided  at  the  expense  of  a  .significant  increase  in  computation 
time  because  of  this  GTD  simulation,  particularly  for  a  large  array  antenna. 


2.3  STP  SIMULATION  EFFICIENCY 

In  the  process  of  enhancing  the  steady-state  STAP  simulation  to  provide  the  new 
capabilities  that  have  been  described,  the  code  has  been  reorganized.  Each  of  the  processing 
options,  configurations  for  each  dimension  of  DOF,  interference  signal  scenarios  (jammers 
and  clutter),  plus  various  aircraft-antenna  interactions  are  described  by  separate  callable 
subroutines.  This  enables  the  simulation  to  be  readily  applied  to  many  radar  scenarios.  It  has 
also  enabled  the  CPU  running  time  for  the  code  to  be  sped  up,  which  can  be  significant  for 
large  arrays.  An  example  of  this  improved  CPU  efficiency  is  shown  in  Figure  17. 
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Number  of  Antenna  Columns 


SECTION  3 


FINITE  SAMPLE  PROCESSING  CAPABILITY 


The  data  received  from  a  pulse-Doppler  radar  is  represented  by  a  non-stationary  train  of 
pulses.  For  the  application  of  STAP,  these  pulses  must  be  recorded  coherently  (e.g.,  as  I  &  Q 
data  for  each  pulse).  An  estimated  covariance  matrix  is  created  by  averaging  the  covariances 
over  a  finite  number  of  space  (adjacent  range  cells)  and  time  (pulses)  samples  from  the 
data  [8].  (This  differs  from  the  steady-state  simulation  where  the  covariance  matrix 
represents  an  average  based  upon  an  infinite  number  of  samples  from  each  clutter  cell.) 
Adapted  weights  are  obtained  by  multiplying  the  inverse  of  this  estimated  covariance  matrix 
by  the  steering  vector.  These  weights  can  now  be  applied  to  the  data  samples  corresponding 
to  the  ranges  and  pulses  over  which  the  weights  were  averaged. 

The  finite  sample  proce.ssing  code  that  was  developed  earlier  [8]  has  been  modified  to 
provide  more  flexible  processing  options.  A  flow  chart  of  the  modified  code  is  shown  in 
Figure  18.  It  begins  with  an  input  consisting  of  a  finite  train  of  I  &  Q  data  samples  that  are 
the  outputs  received  by  the  array  antenna.  These  received  data  samples  can  be  obtained  from 
measured  data  or  from  a  simulation  to  create  the  finite  train  of  received  data  samples. 

A  cube  of  data  is  created  from  this  finite  number  of  data  samples.  The  axes  of  the  cube 
are  the  spatial  DOF  (horizontal  and  vertical  element  outputs  of  the  array),  temporal  DOF 
(pulses),  and  the  contiguous  range  cell  samples.  Options  for  preprocessing  transformations 
of  the  data  are  then  provided:  Spatial  transformations  from  element  space  to  either  beam- 
space  or  subarray  space  and  temporal  transformations  from  pulses  to  Doppler  frequency. 

The  inverse  of  the  estimated  covariance  matrix  obtained  from  all  or  a  subset  of  the  finite 
sample  data  cube  is  multiplied  by  the  appropriate  steering  vector  to  create  the  adaptive 
weights.  These  weights  are  then  applied  to  the  data  cube  to  obtain  the  proce.ssed  output 
.signal.  The  adaptive  processing  options  include  N-pulse  STAP  followed  by  Doppler 
processing  (where  N  is  less  than  the  total  number  of  coherent  pulses  recorded)  and  higher- 
order  Doppler  frequency-factored  processing  [9].  The  available  output  performance 
measures  are  the  interference  signal  level  and  the  SIR  versus  range  or  Doppler  frequency.  To 
illustrate  this  finite  sample  processing  capability,  a  data  file  from  the  Naval  Re.search 
Laboratory  eight-element  airborne  UHF  array  was  proce.ssed.  The  SIR  versus  range 
performance  with  no  processing  and  with  three-pulse  STP  is  shown  in  Figure  19.  For  this 
example,  the  target  signal  at  each  range  is  made  proportional  to  the  steering  vector  and  is  not 
a  function  of  elevation  angle.  In  general,  to  calculate  the  SIR  when  the  target  signal  is  not 
already  in  the  data  (e.g.,  to  insert  a  target  signal)  will  require  the  development  of  a  simulated 
target  signal  as  a  function  of  elevation  gain  pattern  and  Doppler  frequency. 
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Figure  18.  Finite  Sample  Processing  Option 


26 


SIR  Versus  Range  for  Quiescent  Processing  and  Three-Pulse  STP 

for  DLR37 


SECTION  4 


GRAPHICAL  USER  INTERFACE 


A  GUI  has  been  provided  as  an  intuitive  means  for  running  the  software  and  displaying 
its  output  [10],  Interaction  using  the  interface  to  the  software  is  mouse-based  with  pull-down 
menus  describing  output  options.  The  GUI  prompts  the  user  for  an  interactive  editing  and 
verification  of  an  input  file,  runs  the  software,  and  then  offers  several  options  for  viewing  and 
storing  the  output.  This  is  illustrated  by  the  pull  down  menu  in  Figure  20  for  providing  the 
options  for  viewing  the  output  data. 

The  interface  is  written  in  a  relatively  new  interpreted  scripting  language  developed  by 
J.  Ousterhout  at  U.  C.  Berkeley  called  Tool  Command  Language  and  Toolkit  (Tcl  and  Tk). 
This  language  was  chosen  for  several  reasons,  including  its  cross-platform  functionality  and 
the  short  time  and  simplicity  with  which  applications  can  be  developed  and  modified  using  it. 

The  cross-platform  functionality  of  the  GUI  and  the  plotting  software,  GNUPLOT,  (a 
freeware  package  developed  at  Dartmouth)  that  has  been  implemented  with  it  was  important 
in  their  selection  since  the  simulation  runs  on  both  VMS  and  UNIX  platforms.  Tcl  uses  the 
same  interpretation  method  for  both  types  of  operating  systems.  Also,  since  the  script  is 
easily  modifiable  to  accommodate  any  plotting  software  package  that  can  be  loaded  as  a 
series  of  plotting  instructions  at  the  command  line,  other  software  for  creating  plots  can  also 
be  used  for  specific  platforms  for  which  they  are  available,  such  as  MATLAB. 
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SECTION  5 


DEMONSTRATION  OF  STEADY-STATE  SIMULATION 


The  capability  provided  by  the  steady-state  simulation  was  demonstrated  by  its  use  for 
the  Joint  Airborne  Early  Warning  (AEW)  STAP  Requirements  Study.  Simulated  output  I  & 
Q  data  cubes  for  an  airborne  radar  with  a  high  PRF  were  created  from  the  covariance  matrix 
of  the  steady-state  simulation.  A  Cholesky  decomposition  is  applied  to  the  steady-state 
covariance  matrix  for  each  range.  The  lower  triangular  matrix  that  is  created  is  multiplied  by 
a  complex  Gaussian  variance  vector  with  zero  mean  to  obtain  the  I  &  Q  data  samples  for 
each  range.  This  is  repeated  for  each  of  the  unambiguous  ranges.  The  size  of  the  data  cubes 
was  32  elements  x  256  pulses  x  140  ranges  (the  unambiguous  range  inteiwal).  Because  of  the 
design  of  the  antenna,  this  required  a  simulation  with  both  vertical  and  horizontal  spatial 
DOF  in  the  array  aperture.  An  example  of  the  spectrum  obtained  from  two  data  cubes 
created  from  two  Monte  Carlo  runs  for  the  same  scenario  is  shown  in  Figure  21.  The  average 
of  a  large  number  of  these  spectra  approaches  the  .steady-state  spectrum. 

A  second  demonstration  is  ba.sed  upon  a  MCARM  experiment  scenario  described  by  the 
radar  and  environmental  parameters  in  Table  1.  The  clutter-plus-noi.se-to-noise  ratio 
obtained  from  two  independent  simulations  of  the  scenario  are  compared  to  provide  a  basis  of 
confidence  in  the  steady-state  adaptive  simulation  described  in  this  report  [5].  The  other 
simulation,  which  only  provides  performance  prediction  of  an  airborne  radar  without 
adaptive  processing  and  with  idealized  Doppler  filters  (no  “spilling”  between  adjacent  bins), 
has  been  compared  closely  with  measured  data.  One  unambiguous  75-km  range  obtained 
from  the  steady-state  adaptive  simulation  is  shown  in  Figure  22a.  The  performance  for  the 
same  scenario  obtained  from  the  second  simulation  is  shown  in  Figure  22b  for  two  adjacent 
75-km  unambiguous  range  intervals.  In  Figure  22b,  the  antenna  scanned  to  -14.5  degrees  in 
azimuth  and  the  2()()()-Hz  extent  of  the  Doppler  frequency  spectnim  is  represented  by  200 
Doppler  “bins.”  Also,  the  color  scale  of  Figure  22b  plot  is  more  compre.ssed  than  for 
Figure  22a.  The  high  level  of  clutter  signal  calculated  in  the  main  beam,  u.sing  both 
simulations,  is  caused  by  the  high  constant  sigma  model  for  clutter  cho.sen  for  this  scenario. 

The  adaptive  suppression  of  received  clutter  is  demonstrated  with  the  same  scenario  using 
the  steady-state  STAP  simulation.  Figure  23a  is  a  line  plot  of  the  same  unadapted 
performance  shown  in  Figure  22a.  After  two-pul.se  space-time  adaptive  filtering,  the  CNR  is 
significantly  reduced  as  shown  in  the  line  plot  in  Figure  23b. 
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Range  #43  =  151.5331  nmi  =  280.6393  km 
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Figure  21.  Examples  of  Finite  Sample  Clutter  Spectra  Created  From  Simulation 

With  Two  Spatial  Dimensions  of  DOF 


Table  1.  Input  Parameters  for  the  MCARM  Example 


Transmit  Power  (kW)  14.4 

Number  of  Antenna  Columns  12 

Number  of  Antenna  Rows  8 

Frequency  L-Band 

PRF  (Hz)  2()()().0 

Radar  Bandwidth  (MHz)  1.0 

Azimuth  Scan  Angle  (Degrees)  14.5 

Elevation  Scan  Angle  (Degrees)  2.5 

Platform  Altitude  (Feet)  492 1 .5 

Platform  Velocity  (nmi/sec)  219.7 

Platform  Crabbing  None 

Uncompressed  Pulsewidth  (fisec)  lOO.O 

Fractional  Wavelength  Spacing  —  Columns  0.454 

Fractional  Wavelength  Spacing  —  Rows  0.588 

Noise  Figure  (dB)  2.5 

System  Losses  None 

Additional  Losses  on  Target  Only  (dB)  4.5 

Atmospheric  and  Lens  Losses  Included 


Table  1.  Input  Parameters  for  the  MCARM  Example  (Concluded) 


Cancellation  Ratio  (clB) 

Clutter  Model 

Mean  Clutter  Level,  Gamma  (dB) 

Internal  Clutter  Motion 
Number  of  Near-Field  Scatterers 
Number  of  Coherently  Integrated  Outputs 
Number  of  Doppler  Filters 
Doppler  Weighting 
Azimuth  Transmit  Weighting 
Elevation  Transmit  Weighting 
Azimuth  Receive  Weighting 
Elevation  Receive  Weighting 
Number  of  Receive  Elevation  Subairays 


60.0 

Constant  Sigma 
-9.0 
None 
None 
200 
200 

100  dB  Dolph-Chebyshev 
Uniform 
Uniform 

35  dB  Taylor,  N  =  4 
Uniform 
2 


Number  of  Elements  Per  Receive  Elevation  Subarray  4 
Target  Radar  Cross-Section  (m^)  2.0 


Number  of  Non- Adaptive  Pulses  Processed  1 

Number  of  Adaptive  Pulses  Proce,ssed  2 

Steering  Vector  Doppler  Half  Blind  Speed 
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b.  Reference  Non-Adaptive  Simulation 


Figure  22.  MCARM  Scenario,  Demonstration  of  Non-Adaptive  Performance  Using 

Two  Independent  Simulations 
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Non-Adapted  Clutter-Plus-Nolse-to-Noise  Ratio  (dB) 


•III*. 
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Doppler  (Hz) 


a.  Non-Adapted  Received  Clutter 


Adapted  Clutter-Plus-Noise-to-Nolse  Ratio  (dB) 


•iiif 


1000  (Km) 

Doppler  (Hz) 

b.  Received  Clutter  with  2  Pulse  STP 


Figure  23.  MCARM  Scenario,  Demonstration  of  Clutter  Reduction  Usinj 

Steady-State  STP  Simulation 


SECTION  6 


APPLICATIONS  FOR  SIMULATION  AND  PROCESSING  SOFTWARE 


The  capability  and  flexibility  that  have  been  added  to  this  software  will  enable  it  to  be 
usefully  applied  to  many  current  and  planned  programs  for  predicting  performance  and  for 
off-line  processing  of  measured  data.  These  are  outlined  briefly. 

•  MCARM  experiment 

process  measured  data 

evaluate  alternative  processing  architectures 

effect  of  near-field  scattering  for  antenna  platform 

•  predict  performance  (increa.se  of  CNR)  when  beam  is  steered  toward  wing 
proce.ss  measured  clutter  data  with  beam  steered  toward  wing 

•  evaluate  options  for  processing  architectures  to  suppress  near-field  scattering 

•  insert  near-field  scattering  effects  into  clutter  data  measured  using  free  space 
antenna 

-  predict  performance  of  cuirent  and  real-time  MCARM  experiments  and  identify 
beneficial  data  gathering  scenarios 


•  Demonstrate  increased  surveillance  capability  of  slow,  small  targets  for  AW  ACS  radar 
using  STAR 

reconfigured  MCARM  array  (8  rows  x  16  columns,  two  to  four  subarrays/row, 
horizontal  polarization)  to  provide  more  elevation  DOF  and  fit  in  same  radome  on 
BACl-11 

•  model  flight  and  .scan  configurations  to  provide  “J  hook”  characteristics  in 
measured  clutter  data 

•  evaluate  proce.ssing  techniques  ba.sed  upon  row  outputs  for  suppressing  clutter 

evaluate  sidelobe  canceler  architectures  for  AW  ACS  using  auxiliary  radiators 
around  periphery  of  airay  and  guard  channel 
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•  ARPA  mountain  top  data  reduction 

process  measured  data 

inseit  effects  of  near-field  scattering  into  measured  data  and  evaluate  processing 
techniques 

predict  perfonnance  of  future  airborne  experiment  using  RSTER  radar 

•  Rome  Laboratory’s  Bistatic  Adjunct  Surveillance  System  (BASS) 

predict  performance  with  alternative  STAP  architectures 

process  measured  or  simulated  I  &  Q  data  from  receiving  array  outputs 

•  Predict  performance  of  future  airborne  phased  array  radar  candidates  and  evaluate 
alternative  STAP  architectures 
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SECTION  7 


CONCLUSION 


Enhancements  have  been  added  to  our  STAP  steady-state  simulation  and  processing 
software  for  an  airborne  phased  array  radar  that  provide  more  flexible  models  of  the  antenna 
and  processing  architecture  and  more  realistic  models  for  the  simulation  of  the 
electromagnetic  environment.  Also,  a  GUI  has  been  created  that  makes  it  simpler  to  control 
the  input  parameters,  to  examine  input  data  files  and  to  provide  a  flexible  choice  of  output 
performance  parameters  and  output  formats.  As  a  result,  the  simulation  can  be  more 
effectively  used  to  predict  the  performance  of  airborne  radars  and  examine  the  potential 
benefits  of  using  adaptive  STAP  with  these  radars.  Also  measured  or  simulated  finite  sample 
data  can  be  processed  to  determine  the  viability  of  alternative  processing  architectures. 

These  new  features  and  capabilities  have  been  demonstrated  and  potential  applications 
identified. 
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AEW 

ARPA 

AWACS 

BASS 

CNR 

CPU 

DOF 

GTD 

GUI 

JNR 

MCARM 

NECBSC 

PRF 

PRI 

RL 

SIR 

SNR 

STAP 

STP 

Tcl  and  Tk 


GLOSSARY 


Airborne  Early  Warning 
Advanced  Research  Projects  Agency 
Airborne  Warning  and  Control  System 

Bistatic  Adjunct  Surveillance  System 

clutter-plus-noise-to-noise  ratio 
central  processing  unit 

degrees  of  freedom 

geometric  theory  of  diffraction 
graphical  user  interface 

jammer-plus-noise-to-noise  ratio 

Multichannel  Airborne  Radar  Measurements 

Numerical  Electromagnetics  Code-Basic  Scattering  Code 

pulse  repetition  frequency 
pulse  repetition  interval 

Rome  Laboratory 

signal-to-interference  ratio 
signal-to-noise  ratio 
space-time  adaptive  processing 
space-time  processing 

Tool  Command  Language  and  Toolkit 
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MISSION 

OF 

ROME  LABORATORY 


Mission.  The  mission  of  Rome  Laboratory  is  to  advance  the  science  and 
technologies  of  command,  control,  communications  and  intelligence  and  to 
transition  them  into  systems  to  meet  customer  needs.  To  achieve  this, 
Rome  Lab; 


a.  Conducts  vigorous  research,  development  and  test  programs  in  all 
applicable  technologies; 

b.  Transitions  technology  to  current  and  future  systems  to  improve 
operational  capability,  readiness,  and  supportability; 

c.  Provides  a  full  range  of  technical  support  to  Air  Force  Materiel 
Command  product  centers  and  other  Air  Force  organizations; 

d.  Promotes  transfer  of  technology  to  the  private  sector; 

e.  Maintains  leading  edge  technological  expertise  in  the  areas  of 
surveillance,  communications,  command  and  control,  intelligence,  reliability 
science,  electro-magnetic  technology,  photonics,  signal  processing,  and 
computational  science. 


The  thrust  areas  of  technical  competence  include;  Surveillance, 
Communications,  Command  and  Control,  Intelligence,  Signal  Processing, 
Computer  Science  and  Technology,  Electromagnetic  Technology, 
Photonics  and  Reliability  Sciences. 


